function [ sigma, gamma_0, gamma_coll, gamma_pol,gamma_own, gamma_inc, gamma_N, gamma_2N]=unpack_params(params)


global model_type


sigma=params(1,1);
gamma_0=params(2,1);

gamma_coll=0;
gamma_pol=0;
gamma_N=0;
gamma_2N=1;
gamma_own=0;
gamma_inc=0;

if strcmp(model_type,'NoDem')==1 || strcmp(model_type,'NoDemD')==1
        gamma_N=params(3,1)  ;
        gamma_2N=params(4,1)  ;
elseif strcmp(model_type,'CollPerc')==1 || strcmp(model_type,'CollPercD')==1  % college percentage in utiilty
        gamma_coll=params(3,1);
        gamma_N=params(4,1)  ;
        gamma_2N=params(5,1)  ;
elseif strcmp(model_type,'CollPercPol')==1  || strcmp(model_type,'CollPercPolD')==1 || strcmp(model_type,'CollPercPolDOD')==1 % college and percent dem
       gamma_coll=params(3,1)  ; 
       gamma_pol=params(4,1)  ; 
       gamma_N=params(5,1)  ;  
       gamma_2N=params(6,1)  ;
elseif strcmp(model_type,'CollPercPolOwn')==1 || strcmp(model_type,'CollPercPolOwnD')==1 % college and percent republican
       gamma_coll=params(3,1)  ; 
       gamma_pol=params(4,1)  ; 
       gamma_own=params(5,1)  ;        
       gamma_N=params(6,1)  ;  
       gamma_2N=params(7,1)  ;
 elseif strcmp(model_type,'CollPercPolInc')==1 || strcmp(model_type,'CollPercPolIncD')==1 % college and percent republican
       gamma_coll=params(3,1)  ; 
       gamma_pol=params(4,1)  ; 
       gamma_inc=params(5,1)  ;        
       gamma_N=params(6,1)  ;  
       gamma_2N=params(7,1)  ;      
end

%gamma_2N=1;